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SUMMARY 

Within the limitations of linearized potential theory for super- 
sonic flow, solutions have been obtained for the damping in roll of 
triangular, trapezoidal, rectangular, and two swept-back plan forms. 

The results indicate that as the aspect ratio is increased the 
limiting value of the damping in roll for trapezoidal (with a finite, 
fixed rake angle) and rectangular plan forms is equal to the value 
for two-dimensional flow which is twice the limiting value for trian- 
gular plan forms. For the swept— back plan forms having the Mach 
cone and leading edge coincident or very nearly coincident, the 
damping in roll even exceeded the value for two-dimen'sional flow. 

In addition, an investigation of the effect of reversal of the 
plan— form position relative to the stream direction was made for 
the majority of the plan forms considered and the results show 
that this reversal had no effect on the value of the damping— in- 
roll stability derivative. 


INTRODUCTION 

There are a number of methods available for determining 
supersonic— flow load distributions on lifting surfaces by means of 
linearized potential theory. Application of any of these methods 
varies in detail and no individual method can be considered as being 
best suited for use in obtaining solutions for arbitrary Mach-cone 
plan— form configurations. 

The relatively simple problem of determining the induced flow 
field and the surface shape to support an arbitrarily prescribed 
load distribution always can be solved using doublet distributions 
and a surface integration. This procedure provides an explicit 
expression of the doublet sheet potential. (For a detailed descrip- 
tion see reference 1.) On the other hand, if the loading over a 
wing is to be determined from a knowledge of the downwash distri- 
bution required to make the streamlines conform to the shape of 
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the lifting surface, the doublet distribution and surface integra- 
tion method again can be applied directly, provided the wing has 
so-called supersonic leading edges, trailing edges, and tipso 1 When 
subsonic edges or tips are considered, however, the mathematical 
development involves the solution of an integral equation in order 
to get results comparable to those obtained by the relatively simple 
surface integration used for wings with supersonic edges . The 
solution of an integral equation is generally a rather long and 
tedious process and can be circumvented for certain plan forms with 
subsonic leading edges or tips by a method presented in reference 2 0 

In this investigation to determine the linearizedr-potential— 
theory load distribution due to roll, only the surface integral 
methods of references 1 and 2 have been used. Other methods that 
have been applied successfully to obtain angle— of -attack loadings, 
such as the conical flow and the doublet— line methods, are also 
applicable . 

The plaii forms considered in this investigation are the follow^- 
ing: (l) triangular with subsonic leading edges and with supersonic 

leading edges; (2) trapezoidal with all possible combinations of 
raked in, raked out, subsonic or supersonic tips (fig. l) ; (3) 
rectangular; and (b) two swept-back plan forms developed from the 
triangular wings (fig, 2). Moreover, all but the swept— back plan 
forms and the triangular plan form with subsonic leading edges were 
analyzed with the stream direction reversed so that the leading edges 
were interchanged with the trailing edges. 

The load distributions due to rolling obtained for the wings 
investigated were subsequently integrated to determine the rolling 
moment and the damp ing— in-roll stability derivatives. All results 
have been presented in stability— derivative form. 


SYMBOIS AND COEFFICIENTS 
x,y,z Cartesian coordinates 

u perturbation velocity along the positive X-axis 

w perturbation velocity parallel to Z-axis (positive downward) 


1 A supersonic edge is an edge for which the angle of inclination 
from the plane of symmetry is greater than the Mach cone angle. 
The inverse of this relationship defines a subsonic edge. 


free— stream velocity 

span of wing measured normal to plane of symmetry 
root chord of wing 

over-all length of swept-back wing (3ee fig. Z) 
area of wing 


aspect ratio 
density in the free stream 
free-stream dynamic pressure (^pV 2 ) 

pressure differential across wing surface, positive upward 
rolling moment about X-axis 



(i3&) 


rolling-moment coefficient 
rate of roll, radians per second 
damping-in— roll stability derivative 
free— stream Mach number 


(. 

\d(pb/2V) 


) 


J : m 2 -i 

Mach angle 


arc tan 


tip rake angle measured from line parallel to plane of 
symmetry in plane of wing 


tangent of 5 


tan & 
tan h 


complete eliptic integral of the second kind with 
modulus yl— 

complete el liptic integral of the first kind with 
modulus t/l-3 2 


perturbation velocity potential 


k 
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METHOD 


The thin-alrf oil— theory boundary condition for a rolling wing 
is a linear spanwise variation of angle of attack which corresponds 
to a linear spanwise distribution of the vertical induced velocity w. 
The problem, therefore, is to find the load distribution that will 
satisfy this coiidition within the wing boundaries and conform in all 
other respects to a proper solution of the linearized differential 
equation for supersonic flow. A detailed development and discussion 
of the surface integral method for obtaining the load distribution 
for a rolling wing can be found in reference 1. The presentation of 
the method in this report is merely an outline of the operations 
involved in obtaining the load distribution from the known boundary 
conditions. 

The loading at any point is proportional to the perturbation 
velocity u parallel to the longitudinal axis of the wing (the 
X-axis in figs. 1 and 2). The following simple relationships from 
linearized potential theory 


can be used to obtain u and subsequently the load distribution in 

AP 

terms of — once the general expression for the vertical induced 
velocity w required in equation (l) is known. 

To convert the load distribution into rolling moment about the 
X-axis the following integration must be performed: 



w dz 


( 1 ) 


AP _ 4u 

T " v 


( 2 ) 
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The stability derivative coefficient C is then 



8 C 2 2VL 


<j SI)2 p 


The general expression for w must be known in the region 
directly above the wing extending from the wing to the Mach cone 
envelope. The procedures for obtaining this expression are divided 
into the two following categories for discussion. 


An application of Green's theorem, discussed in reference 1 , 
provides the means for determining a general expression for w if 
w is specified in the plane of the wing over the area bounded by 
the envelope of the Mach forecones stemming from the trailing 
edges. The means used amounts to an integration over a doublet 
sheet for which the strength is proportional to the specified values 
of w on the X— Y plane. For the supersonic-edged plan forms, there- 
fore, the loading solution can be obtained readily, since w is 
equal to zero in the bounded area of the X-Y plane except on the 
wing itself where w is equal to the prescribed boundary— condition 
value . 


If the plan form has a subsonic leading edge or tips, the air 
can be disturbed in the region between the foremost Mach cone and 
the wing surface. In this region, therefore, the specification of 
w on the X-Y plane requires in general that an integral equation 
must be solved to obtain the general expression for w. Reference 2, 
however, provides an alternative procedure for certain types of sub- 
sonic leading— edge or tip plan forms whereby the solution of an 
integral equation may be avoided. The solutions for these plan forms 
then are .obtained by a process of integration similar to the 
procedure used for plan forms with supersonic tips and leading edges. 


The stability-derivative results are divided into the trian- 
gular plan form, the trapezoidal and rectangular plan form, and the 


Supersonic Leading Edges and Tips 


Subsonic Leading Edges or Tips 


DISCUSSION OF RESULTS 
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swept-back plan— form categories for discussion. Expressions for the 
stability— derivative coefficients are given. Curves showing the 
variation of these coefficients with aspect ratio are presented in 
figures 3, and 5. Expressions for the loading function AP/q 
are given in Appendix A. 


Triangular Plan Forms 

The expression for C^p for the triangular plan form 2 with 
supersonic leading edges is 



and for the triangular wing with subsonic leading edges it is 


= 


-ite 


k( 2^ e --£1- 
^ 1-02 1-02 


K 


-) 


where E is the complete elliptic integral of the second kind with 
modulus ^1— 0 2 , K is the complete elliptic integral of the first 

kind with modulus Jl-Q 2 , and 0 is a parameter that indicates the 
relative positions of the Mach cone and leading edge. For the sub- 
sonic leading-edge plan forms the value of 0 is less than 1 and 
the leading edge is swept behind the Mach cone. The variation of 
3Cip with the aspect ratio parameter 0A is presented in figure 3 
for both types of triangular plan forms. 


Trapezoidal and Rectangular Plan Forms 

There are two general types of trapezoidal plan forms, defined 
by having the tip either raked in or raked out from the leading 
edge. The tip in either case may be classified as subsonic or super- 
sonic depending on whether the Mach-cone angle is greater or less 
than the rake angle. Thus there are four basic Mach— cone plan— form 


^Triangular plan forms with the point forward will be referred to 
as "triangular plan forms" in contrast to the term "inverted 
triangular plan forms" which will be used in reference to base 
forward triangles. 
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configurations. For three of these configurations the possibility 
of the overlapping of the Mach cones from the tips provides a 
secondary configuration for which it is possible to predict the 
loading distribution and the damping due to roll. All the trape- 
zoidal plan forms considered are shown in figure 1. 

The expression for pCj for the trapezoidal plan forms with 
supersonic tips (0>l) P 


was found to apply regardless of whether the tip was raked in, as 
on plan form a, or raked out, as on plan form b. Since for a 
given 6 and a given span the plan forms a and b are effectively 
the same plan form with the stream direction reversed, it is evident 
that Clp is independent of the direction of the stream, provided 
the stream remains parallel to the plane of symmetry of the plan 
form. The same conclusion applies to the subsonic— tipped trapezoidal 
and rectangular plan forms (0<l) and the expression for 0Ci p for 
this type is 


This property of reversibility appears to be rather remarkable 
in view of the fact that the load distribution on a given plan form 
is changed markedly by reversing the stream direction. This result 
also is a direct analogue to von Karmen' s independence theorem 
(reference 3) for drag and, although the drag theorem was developed 
analytically to apply to all plan forms, the independence of the 
lift and roll characteristics as yet merely has been indicated by 
calculation on a limited number of plan forms. 

It is obvious that the trapezoidal plan form can be reduced in 
span till either a triangular or an inverted triangular plan form 
is obtained. Thus it was possible to check on the reversibility 
property of the supersonic— leading-edg9 triangles quite readily. 

It was not possible in the present investigation, however, to reduce 
the subsonic— tipped trapezoids to triangles corresponding to 
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inverted triangular wings in order to obtain a check on the inde 
pendence theorem for these plan forms. As PA becomes less than 
(1 + 0 ) 2 the Mach cone from one tip is reflected off the other tip 
and the load distribution behind these reflected Mach lines is not 
readily determined. From figure 3 it is apparent that, for values 
of £A less than (l + 0) 2 , the variation of £Ci p with PA 

must reduce in slope as the trapezoid reduces to a triangle 

(PA = 40) in order to yield the value of 3Cz p corresponding to 

the value of pCz p for the triangular plan form. 

It can be seen in figure 3 that as the aspect ratio increases 
the values of pCz p for the trapezoidal and rectangular plan forms 
approach the value for two-dimensional flow - ^ ; whereas the 
triangular plan forms reach a value of — ^ when the aspect ratio 

parameter PA becomes equal to or greater than 4. By letting the 
tangent of the rake angle approach infinity as the span approaches 
infinity, the supersonic— tipped trapezoidal plan forms can be made 
to approach the configuration of an inverted triangular plan form. 
Since for a triangular plan form 

A = 4m 

the limiting value of pC} p for this case becomes (using the expres- 
sion for the supersonic— tipped trapezoidal plan forms) 



and thus checks with the value for the triangular plan form with 
supersonic edges. 


Swept— Back Plan Forms 

The two swept— back plan forms shown in figure 2 are easily 
developed from the triangular wings of supersonic and subsonic lead- 
ing edges by cutting out a small triangular area from the rear of the 
triangular plan forms. Other swept— back plan— form configurations 
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are also amenable to calculation but only the results for the two 
plan forms selected are included in this report for purposes of 
comparison. The expression for the damping in roll of the swept-back 
plan form with subsonic leading edges is 


Pda-Ig) 2V = P(lq— L 2 )2V 

qSb pb ^ + 2c r) ~ * 2 ~°r 2 ]]qpb 2 


where Li is the rolling moment for the triangular wing aa a whole 
and L 2 is the rolling moment of the small triangle cut out to 
provide a swept— back plan form. The expressions for Lx and L 2 
are given in Appendix B. 

For the swept— back plan form with supersonic leading edges the 
damping— in— roll coefficient is 

p(L3-L 4 ) 27 3(Ly-L 4 )2V 

P qSb pb m + 2c r ) - l 2 -c r 2 qpb 2 

where L 3 is the rolling moment for the triangular wing as a whole 
and L 4 is the rolling moment of the small triangle cut out to 
provide a swept— back plan form. The expressions for L 3 and L 4 
are given in Appendix B. 

From figure 5 it can be seen that both of these swept>~back 
plan forms provided more damping in roll than their related trian- 
gular plan forms of the same span and Mach— cone leading— edge config- 
uration. This result was anticipated, since the load distributions 
for the triangular plan forms in roll reveal that the magnitude of 
the load increases as the leading edges are approached. At values of 

9 in the immediate vicinity of 1, the swept-back plan forms provide 
greater damping in roll than the trapezoidal or rectangular plan forms 
for comparable values of pA. Moreover, even the value of pCjp 

for two-dimensional flow — j is exceeded in magnitude at PA‘s of 

10 or higher by the swept-back plan forms of the type shown on which 
the Mach cone and leading edge are very nearly coincident (9 very 
nearly equal to l) . 
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CONCLUDING REMARKS 


The results of this investigation based on linearized potential 
theory indicate that for values of PA greater than 4 the trapezoidal 
and rectangular plan forms provide values of -PC^p greater than 

the limiting value j for the triangular plan forms. The limiting 


value of — PCZp for the rectangular and trapezoidal (with a finite, 

fixed rake angle) plan forms is ^ which corresponds to the value 
for a cons tan i>-chord plan form of' 5 infinite span. The value of 

^ for — pCip is also exceeded by some of the swept-^back plan forms. 


The swept-back plan forms investigated, moreover, showed that for 
configurations on which the Mach cone and the leading edge are very 
nearly coincident the value of — PCi for two-dimensional flow 

p * 

%■ is exceeded. 


All but the swept— back plan forms and the triangular plan form 
with subsonic leading edges were investigated with the leading edges 
and trailing edges reversed relative to the stream direction. For 
a given plan form at a given Mach number it was found that this 
reversal, which changed completely the distribution of the load, 
had no effect on the value of PC ip. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX A 

The expressions for the load distribution over the triangular plan 
f orms are as follows : 


e >1 (Supersonic leading edges) Region within the Mach cone 


£P = 4p 

q V 


I 


m(0 2 y-mx) 

*(e 2 -i) 3/2 J 


sin 


— l 


0 2 y-mx 

-e(mx-y)- 


m(02y-tmx) 

n(02_l) 3 /2 


sin 1 


02y+mx 

_ 0(mx+y) 


(02 — l ) 3 / 2 
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Region between Mach cone and leading edge 



0<1 (Subsonic leading edges) 
AP 4p m2yx 



The expressions for the load distributions for the trapezoidal 
plan forms shown in figure 1 and for the rectangular plan form are 
as follows: 

e >i 

Plan form a, PA>4 

AP _ 4 py 

q PV 

Plan form b, pA >4(0+l) 2 /(0+2) 


Region 1 

AP _ 4 py 

q ‘ PV 

Region 2 



sin -1 

L e(xm -y + | -c r ) 



q V it (S 2 -!) 3/2 



my 

— sin -1 
*0 


•1 £ + — 


'e(y-| +mc r) 1 


2(02-1) 3/2 n0 


mx 


20 



n 0 2 (02-l) 
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Plan form c. 
Region 1 
Region 2 
Region 3 


4<3A<4(e+l)2/(0+2) 

same as for region 1 on plan form b 

same as for region 2 on plan form b 

^ same as for region 3 on plan form b 


Region 4 


AP /'ap"\ 

\ -(*£) 

_ ii-fi 

m(-0 2 y-nix- ^ +mc r ) 

sin -1 

0 2 (-y- ^ +mc r )-mx 

q Vq J 

2 V q/ 

1 T l 1 

L *(0 2 -l) 3 / 2 J 

L 0(xm-y+ | -c r ) J 


m (-0 2 y-nix— J + mc r ) my 

2(02-1)3/2 + Se sin 




+mc. 


mx 


my 

20 


m 


Jm 2 x 2 -6 2 (-y— +mc r ) g 

ne 2 (e 2 -i) 


I 


I 


e <i 


Plan form d, 3A > 4/( 2-0 ) 
Region 1 
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Region 2 


AP [ y , 2p(y— fe)+x+0x 

— = — i — cos -1 

q PV L n x(l-e) 


2 /-( Py-g|- +x ) ( py-pjj- +6>x ) 
n(3(l-0) 


Plan form e. 
Region 1 
Region 2 
Region 3 


(1+0)2 <PA <4/(2-0) 

^2 same as for region 1 on plan fonn d. 
same as for region 2 on plan form d. 


AP _ /AP^ 

q ~ \<1 72 




fy [ 2p(-y-^)+x+0x 

u L 


g^Py+P^ -x)(0x-Py-P^) \ 
+ it (3( 1—^ ) J 
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Plan form f, PA >4(0+l) g /(0+2) 


Region 1 

AP 4py 
q PV 


Region 2 


AP 4 p | y _ x f" 2(Py-P^ +0c r )+x-0x ~] 

q PV In 008 L x(l+6) J 


2 ,/(x+Py-P£ +0c r )(ex-py+p| -ec r ) 
n p( 1+0 ) 42 


! 6e g (Py+x)-2e( Py+x+40c r -4p|~) /(x+py-P& +0c r ) \ 

^ 3^P(l+P) 2 * v (0x— Py+P^ — 6c j.) ^ 


Plan form g, (l+0) 2 <PA <4(l+0) 2 /(0+2) 


Region 1 ^ same as for region 1 on plan form f. 


Region 2 ^ same as for region 2 on plan form f. 
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Region 3 


AP f ap'N f LS\ 4 p fy _ x f 2( -gy-f7 + 0 c r )+x- 0 x 

i~ " Vq. /2 A A + C ° 8 i+0 ) J 


2 /(x-Py-p| +0c r )(0x+py+p£ -0c r ) 
*p(l+0) 


60 2 ( Py-x ) +20 ( x-Py+40c r -kp 

/x-Py-P^ +0c r 1 

3np(l+0) 2 J 

Y 01 +Py+P^- -0c r J 


0 = 0 (rectangular) PA >1 


same as for trapezoidal plan form e with 0=0. 

<1 


APPENDIX B 

The expressions for the rolling moments Li, L2, L3, and L4 
for the swept— back plan forms are: 


Li = 


-pqb 4 n 


3 2 V E 


1-02 


K 


i6 
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-8pqm 2 / 

c r 4 e^46 a +i) " 


L 8p 2 ( i-e 2 ) 7 / 2 J 



_ am" 1 

CjjB 


*/Z 2 8 2 — ( Z— c p) 2 
2^3 2 (i-e 2 )3 


31 ! ^-5l 2 c r +lc r +c r 3 -e 2 (9l 3 -13Z 2 c r +72c r 2 +10c r 3 ) 


+ e 4 (9z 3 -nz 2 c r +6zc r 2 -6c r 3 )-3e®( z 3 -z 2 c r ) 


c r 4 e 3 (2e 2 +i3) 
243^ 1-e 2 ) 3 + 



l 3 


-^pqZ^nr 4 

3V0 
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— 2pqm 4 
nV(6-l)3/e 


3Q g +2 ) 

_12(e*-l) ; 


(0 2 -g) 

12 


sin 1 



( 3 e 2 + 2 )(e 2 -k) 

36(6 £-1)3/2 


r (e g + i)(z^ r )3 + (z-c r ) 

L 3c r 3 2c r 2 


f + 

24 J 


sin” 1 

0( 2Z-0 


(9 g -l)(Z-^ r ) 3 t ( Z— c r ) 2 (9 g -2) (39 2 +2) 1 _ x 9 g ( Z-C r )-Z 

3c r 3 + 2c r 2 24 12(9 2 -l) 2 J S n 9c r 


+ 


/z g -9 2 (z-< r ) 2 r_ (z-c r ) g /e 2 -i ^ 

Cj. L 93j. 2 


(Z-c r )(39 2 n-2) 

36c r 


(39 2 +2)(e 2 -4)' 
36(92-1 ) 3 /^ . 


+ 


z 4 


0 2 (Z-* r ) 3 (Z-c r ) 2 

JT 3 2Z* 


(29 g — 3) 

6 


sin' 


0^( Z— Cj»)+Z 
e(2 _ z-c r ) 


+ 


9 g (2-c r ) a _ ( Z— c r ) 2 
3Z 3 2Z 2 


(29 2 -3) 

6 


sin 


0^( Z— c r ) — Z 

0c r 


n(0 2 -l) 3 / 2 ( Z— c r ) y9 S -l /z 2 -9 2 (Z-^ r ) ; 


39 


3Z‘ 


2(0 2 _l)3/2 

30 


cos- 1 > 
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Figure /. - Trapezoidal pfanforms and Mach cone con- 
figurations investigated . 
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Figure 2.- Swept-back p/anforms and Mach cone 
con figuration s investigated. 
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